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It is shown that the two-gap approximation is applicable for describing the dV/dI{V) spectra of 
LuNi2B2C-Ag point contacts in a wide interval of temperatures. The values and the temperature 
dependences of the large and the small gaps in the ab plane and in the c direction were estimated 
using the generalized BTK model^" and the equations of.^^ In the BCS extrapolation the critical 
temperature of the small gap is 10 K in the ab plane and 14.5 K in the c direction. The absolute 
values of the gaps are Ag' = 2.16 meV and Ag = 1.94 meV . For the large gaps the critical 
temperature coincides with the bulk Tc, tJ'""' = 16.8 K, and their absolute values are very close, 
being about 3 meV in both orientations. In the c direction the contributions to the conductivity 
from the small and the large gaps remain practically identical up to 10 11 K. In the ab plane the 
contribution from the small gap is much smaller and decreases rapidly as a temperature rises. 

PACS numbers: 63.20.Kr, 72.10.Di, 73.40.Jn 



I. INTRODUCTION 

Compounds i?eNi2B2C Re is a rare earth metal) have 
a tetrahedral crystal structure (e.g., see^ Fig. 6) with the 
spatial group lA/mrnm^ and the characteristic parame- 
ters: two lattice constants a and c {a — 3.46 A, c = 10.63 
A for LuNi2B2C) and the distance z between the B atom 
and the RC plane. The anisotropy in these compounds is 
determined by the ratio c/a and is rather high (c/a ~ 3).^ 

According to electron structure calculation, i?eNi2B2C 
compounds are three-dimensional metals in which all 
atoms contribute to the density of states at the Fermi 
surfacej^i^ In LuNi2B2C three energy bands crossing 
the Fermi level subdivide the Fermi surface into three 
distinct regions: (1) two spheroidal sheets at F-points 
that are extend along the c axis; (2) two flat P-centered 
square regions with the sides parallel to (100) and (010), 
and finally (3) the largest basic X-centered part mostly 
cylindrical in shape (closely resembling practically two- 
dimensional crystal structure) whose axis is parallel to 
the c axisjSi^ In this sheet there are flat regions at the 
points 0.56±(27r/a), the nesting, typical for the com- 
pounds of this familyi^ It is precisely these features 
that determine the Kohn anomaly near the wave vec- 
tor Qm ~ (0.5, 0, 0)^ and the incommensurate ordering 
with Qm ~ (0.5,0,0) in AFM compounds. The part 
of the Fermi surface responsible for nesting makes only 
(4.4 ± 0.5%) and shows up a sHghtly increased resistance 
for the current along (100).^- 

As was shown for point contacts with a ballistic transit 
of electrons in a many-band superconductor—, in any di- 
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rection the contribution to conductivity from each band 
is proportional to the area of the Fermi surface projection 
from the corresponding band onto the interface. Since 
the fraction of the Fermi surface responsible for nesting 
is small, the nesting-related features are not obvious in 
the point-contact spectra. 

Thus, the Fermi surface of nickel borocarbides con- 
sists of several bands: it is anisotropic and correspond- 
ing Fermi velocity vp varies within wide limits. In this 
case Tc is determined not by the total density of states 
N{Ep), but by the narrow peak of the density of states 
at the Fermi level (see^. Fig. 12). The density-of-states 
peak is induced by the slow electrons from the flat (nest- 
ing) regions of the Fermi surface with the wave vectors 
(0.5-0.6,0,0).! 

In the normal state the resistivity p is isotropic be- 
cause it is connected with the groups of electrons that 
have rather high velocities vp and do not belong to the 
flat regions of the Fermi surface. The temperature de- 
pendence of resistance p in single crystal YNi2B2C and 
LuNi2B2C also exhibits this practically isotropic metallic 
behavior (seei. Fig. 11). 

The scanning tunnel microscope measurements on 
LuNi2B2C give a rather low value of superconducting en- 
ergy gap along the c axis (4.2 K), ~ 2.2 meV {2A/kTc — 
3.2)MiS. In the ab plane (T ~ 0.3 0.5 K) for YNi2B2C 
the gap data are greatly scattered, from 2.3 meV^ to 
3.5 meViii {2A/kTc = 3.44 H- 5.24). The high ratio inii 
can be attributed to the estimation procedure: the gap 
was estimated from the positions of the maxima in the 
differential conductivity curve. The shape of the curve, 
however, suggests a considerable variation of the gap over 
the Fermi surface, and the cited value of the gap is closer 
to its maximum rather than to its average magnitude. 
For a number of nickel borocarbides the ratio obtained 
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from point contact measurements in the ab planeiati^ is 
2A/kTc = 3.7 3.8. 

In the recently discovered superconductor MgB2 the 
Fermi surface is formed by two groups of energy bands 
having different dimensions, which determines a two-gap 
conductivity^^. Similarly, we can expect a many-gap con- 
ductivity in nickel borocarbide compounds. Note that 
the dHvA data^^-^^ also point to the presence of at least 
two or even three groups of the Fermi velocities evidenc- 
ing in favor of a two or three-band model. Inii the be- 
havior of Hc2 (T) is considered within a two-band model. 
One band has low Fermi velocities, high Hc2{0), T^, A 
(electron-phonon interaction) parameters. However, all 
these characteristics are suppressed in the bound system 
by another band with high Fermi velocities and smaller 
superconducting parameters. In this case the two-gap 
approximation commonly used for two-band supercon- 
ductorsiSi seems quite reasonable. 

Here we report the results on energy gap anisotropy 
that were obtained on the nickel borocarbide supercon- 
ductor LuNi2B2C for the main crystallographic orienta- 
tions - in the ab plane and along c axis. This study is a 
logical extension ofi^. It was shown that the one-gap ap- 
proximation could not provide an accurate description of 
the experimental curves in the low temperatures region 
and use of the two-gap approach was contemplated.^^ 
The largest and the smallest gaps were estimated in the 
ab plane and along c direction at T ^ Tc. Besides, using 
the generalized BTKSU and Beloborodko's^ models, the 
lowest values of the largest gap were estimated in these 
directions as a function of temperatureiiS Owing to the 
high quality of the experimental curves (no " humps" , or 
broad maxima at \eV\ > A), the estimation was possible 
in the "wings", i.e. in the regions of the experimental 
curves where the biases are higher than that at the dif- 
ferential resistance minima {\eV\ « A). 

Below we show the results of the two-band calcula- 
tion for LuNi2B2C obtained within two models^^-^^ for 
experimental curves obtained before^ In our opinion, 
the calculation most closely corresponds to the many- 
band character of the Fermi surface in these compounds. 
The anisotropy is estimated qualitatively in each band 
using the parameter F in the BTK modelpS or 7 in Be- 
loborodko's model. The calculation is compared with 
the results of one-gap fittingiiS 



II. EXPERIMENTAL TECHNIQUE 

The point-contact measurement was performed on sin- 
gle crystal LuNi2B2C (T^ ~ 16.9 K) grown by Canfield 
and Budko using a flux method. The crystals were thin 
(0.1 ~ 0.2 mm) plates with the c axis perpendicular to 
the plane of the plate. The single crystal surface always 
has quite a thick layer on it in which superconductiv- 
ity is either absent or strongly suppressed. Point con- 
tacts between natural crystal faces therefore suffer from 
this disadvantage. For measurement in the ab plane, the 



crystal is usually cleaved, and the point contact is made 
between a metallic counterelectrode and the cleaved sur- 
face. The cleavage perpendicular to the c direction is a 
technical challenge. To do this, the crystal surface was 
cleaned with a 10 15% HNO3 solution in ethanol for 
several minutes. The properly etched crystal has metallic 
luster and its surface is free from colored film. The crys- 
tal was washed thoroughly in pure ethanol and than in- 
stalled in contact-making device. As found subsequently, 
the etching conditions are very important for producing 
high-quality point-contacts.^'^ In the ab plane the mea- 
surements results seem to be indifferent to cleavage or 
etching - they are practically similar in both cases. But 
a deeper analysis (see below) reveals certain distinctions 
caused by defects in the cleavage surface. The other elec- 
trode was made of high-purity silver. The point contacts 
in the ab plane were fabricated between the edge of the 
prism-shaped Ag electrode and freshly cleaved (etched) 
face of the single crystal (shear technique^'*). The devia- 
tion from the direction perpendicular to the c axis varied 
within 5 10°. To made a contact in the c direction, 
the "needle-anvil" geometry was used. The needle ra- 
dius was 1-3 microns. The temperature was measured 
with a special insert (similar to that inS^). 

The point contact resistance varied typically from sev- 
eral Ohms to tens of Ohms. To attain more detailed data, 
we selected point-contacts with the greatest possible tun- 
neling which was judged from the presence of a differen- 
tial resistance maximum at zero bias and the highest non- 
linearity of the I-V curves at biases ±A typical for a non- 
perturbed superconducting surface in the contact neigh- 
borhood. Unfortunately, we were able to take a com- 
plete set of curves in the whole range Tmin ~ 1.5K — > Tc 
only on a few contacts. Because of high resistance and 
long-term (over 10-12 hours) measurements of temper- 
ature, many of contacts were broken down. The mea- 
surements were made with the nearly equal temperature 
steps. There were taken 47 curves in the c direction and 
41 curves in the ab plane. 



III. DATA PROCESSING 

Some curves of the temperature series for the first 
derivatives dV/dl of the LuNi2B2C-Ag point contacts in 
the c direction and in the ab plane are shown in Fig.^ 
The contact diameters were calculated by Sharvin's equa- 
tion using iJg'' = 22.5rj, i?g = 4517, pi = 11.25-10-i617to2 
of^"^ and allowing for the presence of a potential barrier 
in the constriction region: Zat — 0.7, Zc — 0.55. The 
obtained diameters dab — 13.7nm and dc — 8.5nm are 
very close in the order of magnitude to the coherence 
length in this compound (^ — Q.5nmM>). The measured 
curves were symmetrized dV/dlsym = l/2[dV/dI{V) + 
dV/dI{—V)] and normalized to the normal state at T > 
Tc. 

Two approaches were used to compare the theoretical 
and experimental results obtained in the one- and two- 
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Here 7=1/ (T5A) is the pair breaking parameter, is 
the mean free time during spin-flip scattering at impuri- 
ties. When magnetic impurities are absent, tends to 
infinity and the equations describing the I — V charac- 
teristics in^^ and in^^ coincide. 

The other approach was based on the generahzed 
Blonder-Tinkham-Klapwijk (BTK) model commonly 
used to describe S — c — N point contacts. The model 
allows for finite lifetime of quasiparticles t = h/T deter- 
mined by inelastic scattering, which leads to the broad- 
ening of the density of states in the superconductor—. 
Formally, according to the previous theory, the BTK- 
based results are obtained under condition of strong pair 
breaking (|u| ^ Therefore, in the strict sense, the 

generalized BTK model contains no gap. For any in- 
finitesimal broadening parameter F, at T = the density 
of states near Fermi surface is nonzero. In theorjiSi, the 
order parameter is qualitatively analogous to the pseu- 
dogap in the generalized BTK model. 

In the two-gap approximation, the general conductiv- 
ity of the point contact is taken as a superposition of the 
conductivities from two regions of the Fermi surface hav- 
ing their own gaps. In this case the experimental curves 
were fitted using the expression 
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FIG. 1: Differential resistances of LuNi2B2C-Ag point con- 
tact at different temperatures. To avoid crowding, only sev- 
eral curves are shown. 



gap approximation. First, a model was applied, which 
describes the electrical conductivity of pure S — c — N 
point contacts in the presence of an arbitrarily transpar- 
ent potential barrier at the boundary between the metals. 
The model allows for the finite lifetime of Cooper pairsjSi 
The equations describing the I —V characteristics of the 
point contact within this model are given in^SiSi The 
model, however, interprets the superconducting order pa- 
rameter and the energy gap as distinct magnitudes. The 
BCS order parameter A was found from 



A = A/&tanh(^)/Ze-^ 



(1) 



K is the coefficient, characterizing the contribution 
from the part of the Fermi surface with the gap Ai, S 
is the scaling factor describing the intensity of the ex- 
perimental curve. It is used to make the amplitudes of 
the theoretical and experimental curves equal. The best 
agreement between the shapes of the theoretical and ex- 
perimental curves corresponding to the minimum rms 
deviation in the curve -F'(A) was taken as the main cri- 
terion of fitting in both the models at low and moder- 
ate temperatures. When the temperature increases, the 
curves i^(A) become flatter, which entails much uncer- 
tainly in estimation of A. In this case the scale-factor S 
is of primary importance. It is, as a rule, independent 
of temperature and at low temperatures it can be calcu- 
lated quite accurately through averaging. The details of 
the calculation in the one-gap and two-gap approxima- 
tions are offered in the Appendix. 



IV. CALCULATION 



where u is obtainable by solving the equation 
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The energy gap Aq and the order parameter A are related 



Ao = A (1-7^/3 



3/2 



(3) 



A. One-gap approximation 

Although LuNi2B2C exhibits a two-band kind of super- 
conductivity, its experimental curves can be described in 
the zero approximation with an averaged gap possessing 
the broadening parameter F or the pair-breaking parame- 
ter 7. The temperature dependences of the order param- 
eters (the model o^) and the gaps (BTK calculationSS) 
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FIG. 2: Temperature dependences of the order parameter 
(energy gap) A calculated in the one-gap approximation— 
(open circles) and within the generalized BTK model^" (solid 
circles) for point contact LuNi2B2C-Ag in the ah plane and 
in the c direction. Solid lines show the BCS approximation. 
Here and in all other figures the following symbols are used: 
circles - one-gap or average gap calculation; squares - the 
small gap; triangles - the large gap; closed symbols - BTK 
calculation Ref.i2Q: open symbols - calculation by formulae of 
Ref.ES 



in the ah plane and in the c direction are shown in Fig. 
121 along with the BCS curves. 

In both cases we observe a deviation from the BCS 
curve. In the c direction the deviation is towards higher 
values when the calculation is based on^^ and towards 
lower values in the BTK approximation. As for the ah 
plane, it should be noted that in the low temperature re- 
gion the resistance of the contact varied during the mea- 
surement and became stable only at 4 K. This may be 
the reason why in both models the deviations from the 
BCS dependence in the low temperature region descend 
to lower values. The results of the both approximations 
start to coincide near 8 X in the ah plane and near W K 
in the c direction. Precisely at these temperatures (see 
Fig. 3), the parameters 7 (pair breaking) or F (broad- 
ening) turn to zero in these orientations. 

The ratio F/A or the 7 magnitude can be used for 
semiquantitative estimation of the contribution of the 
two-band structure to the point contact conductivity. 
Since in the c direction the ratio F/A is 1.5-^2 times 
higher, we can expect a stronger two-band effect in it. 
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FIG. 3: Temperature dependences of broadening F (closed 
circles) and pair-breaking 7 (open circles) parameters for 
LuNi2B2C-Ag point contact in the ah plane and in the c di- 
rection in the one-gap approximation. 



Besides, we can conclude that the experimental curves 
demonstrate one-gap superconductivity above 8 in 
the ah plane and above 12 K in the c direction. True, 
the more accurate two-band calculation (see below) gives 
somewhat different results; nevertheless, this technique is 
quite good for semiquantative calculation. Note that in 
the strict sense, neither F nor 7 are meant for describing 
distribution of gaps over the Fermi surface. They allow 
for the finite lifetime of carriers, which only simulates the 
energy gap anisotropy. Their use in this particular case 
is, however, quite justified as it enables an approximate 
modeling of the I —V characteristic in the N — c — S 
point contact. 

In ah and c orientations the gap is anisotropic, though 
observed anisotropy is rather low (see Figs . 1112) 1. Statis- 
tically (several tens of contacts were examined for each 
direction), the positions of the dV/dl-mmima character- 
izing the gap value are 15 20% shifted towards higher 
energies in the ab plane. 

In addition to the deviation from the BCS A(r)-curve, 
the one-gap model also fails to provide a good fitting of 
the shapes of theoretical and experimental dV/dI{V)- 
curves. The one-gap fitting (BTK model^'^) of experi- 
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FIG. 4: Fitting of experimental curves in the one-gap ap- 
proximation (BTK model^") for LuNi2B2C-Ag point contact. 
A few experimental points are shown to avoid overloading. 



FIG. 5: Temperature dependences of the rms deviation char- 
acterizing shape discrepancies between experimental and the- 
oretical (calculation by Ref. [^curves (open circles) and in the 
BTK model Ref.|20| (closed circles)) for LuNi2B2C-Ag point 
contact in the ab plane (a) and c direction (b) in the one-gap 
approximation. 



mental curves is shown in Figs.^] and [S] illustrates the 
temperature dependences of the rms deviations between 
the shapes of the theoretical and experimental curves. 
The best agreement is observed above 7 ^ 8 K. 

We attribute this rather high value to the two-band 
character of the superconductivity. This indicates that 
experimental curves must be fitted in the two-gap ap- 
proximation, as in the case of MgB2 andiS. 

B. Two-gap approximation 

The two-gap fitting of experimental curves in the low 
temperature region is illustrated in Fig. |2J5. The agree- 
ment is seen to be much better than with the one-gap 
approximation (compare with Fig. 

The temperature dependences of the order parameters 
calculated from;21 (open symbols) and the gaps calculated 
within the generalized BTK modcP" (solid symbols) in 
the ab plane and along the c axis are shown in Fig.[7| 
along with the corresponding BCS dependences (solid 
lines). In addition to the small Ai (squares) and large 
A2 (triangles) gaps, Fig. [7| illustrates for a comparison 
an averaged gap Am (circles) which allows for the partial 



contributions to the conductivity of the contact made by 
the Fermi surface areas with different gaps (see formula 
in Fig. El. 

It is interesting to compare the BCS approximation 
of the critical temperatures of the small gap in the ab 
plane and in the c direction: T^''(Ai) = 10 K and 
r^(A2) ~ 14.6 K. The absolute values of the small gaps 
are close (Ag'' = 2.16 mV, Ag = 1.94 mV), the differ- 
ence being no more than 10%. The critical temperatures 
of the large gaps are similar and coincide with the bulk 
Tc = 16.8 K. The values of the large gaps are about 
A2 = 3 mV. Note that both the models give practically 
identical results for the large and the small gaps (open 
and closed triangles and squares). As for the average 
gaps (circles), their values in the low temperature region 
are dependent on the model applied. This is because of 
different correlations between the partial contributions 
from the large and the small gaps (see Fig. 

It thus turns out that the higher F (or 7) in the c direc- 
tion (where the two-band character is more pronounced) 
obtained in one-gap approximation is determined not by 
the difference between the values of the large and the 
small gaps but by the difference between the partial con- 
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FIG. 6: Two-gap approximation (BTK model Ref.|23 of ex- 
perimental curves for LuNi2B2C-Ag point contacts. To avoid 
crowding, only a part of experimental points is shown. 

tributions from these gaps. In the ab plane the contribu- 
tion of the small gap (open squares, Fig. |S| upper panel) 
is smaller than in the c direction and it decreases rapidly 
as the temperature rises. In the c direction the contribu- 
tions from both gaps are very close (in average between 
the two modelsSSiSi is about 50% for each) up to 10-;- 11 
K. (Note that in the BTK calculation^*^ the contribu- 
tion from Ai is over 50%; in the calculation by model 
o^ the Ai contribution is below 50% and in average 
the contribution from Ai is close to that from A2 at low 
temperatures.) 

The temperature dependence of 7 (pair breaking), or F 
(broadening of quasiparticle levels) obtained in the two- 
gap approximation are shown in Fig. |51 

Note that the parameters F and 7 are considerably 
higher for the small gap both in the ab plane and in the 
c direction. Their ratios to the corresponding gaps are 
larger too: (Fi/Ai) = 0.185 > (F2/A2) = 0.021 in the 
ab plane and (Fi/Ai) = 0.333 > (F2/A2) = 0.1 in the 
c direction. As stated above, we consider that these ra- 
tios simulate qualitatively to what extent the gap is dis- 
tributed over the corresponding sheet of the Fermi sur- 
face. We can thus assume that at low temperatures the 
value of the gap changes very little for the part of the 
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FIG. 7: Temperature dependences of the order parameter 
(energy gap) A calculated in the two-gap approximation us- 
ing the equation of Ref.|^ (open circles) and the BTK model 
Ref.|23 (solid circles) for LuNi2B2C-Ag point contacts in the 
ab plane (a) and c direction (b). Solid lines - BCS extrapo- 
lation. In addition to the large A2 (triangles) and small Ai 
(squares) gaps, the figure includes medium gaps Am (circles) 
that were obtained taking into account the partial contribu- 
tion to the contact conductivity from the Fermi surface re- 
gions with the large and small gaps (see the formula in the 
figure and Eq.(4) in the text). 

Fermi surface with the larger gap (A2). In the small- 
gap region (Ai), the gap is spread over a considerably 
broader range of energies at low temperatures. Thus 
the two-gap approximation, even being more accurate, 
cannot describe the experimental curves in the whole T- 
intcrval. This is evident (see Fig. 10) in the temperature 
dependences of the rms deviations F between the shapes 
of the fitted and experimental curves (even for two-gap 
calculation). 

Although the errors observed (the F magnitude) are 
much lower (about half as high at the lowest tempera- 
ture) than in the one-gap case (Fig. [SJ, they increase 
appreciably below T — 4.5 5 K. This may occur be- 
cause the three-band conductivity comes into play, which 
must not be ruled out completely for this compound. 
Proceeding from the significantly higher ratio Fi/Ai at 
low temperatures, we can treat the small gap in its turn 
as a superposition of two different gaps with different 
Tc- These gaps are close to each other and their energies 
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FIG. 8: Temperature dependences of the partial contribution 
to the contact conductivity from the small gap (calculation 
by the formula of Ref.ll^ and within the BTK model Rei^ 
Also, see Eq.4 in the text. 

overlap. The temperature dependence of the error (F) 
suggests that we can expect Tc ^ 5 K for the smallest 
gap in the BCS approximation. This is the lowest es- 
timate. However, because r2 becomes zero at 6 K (a 
plane) and 8 K {c direction) (Fig. |3, it is possible that 
precisely these temperatures are closest to Tc of the gap 
in the third band. 



V. DISCUSSION 

In our experiments wc have not succeeded in detecting 
noticeable anisotropy of the energy gap in the ab plane 
and in the c direction. It was no more than 20% in the 
average-gap approximation. The reason may be a rough 
cleavage surface, which made it difficult to fix the contact 
axis precisely along a direction where the gap minimum 
is predicted. Therefore, following the common practice in 
point contact spectroscopy, we tried to select point con- 
tacts with the highest superconducting parameters (gaps, 
I — V nonlinearities). Since in the ab plane the maxi- 
mum gap is predicted in the (110) direction, our contacts 
probably were made to correspond most closely to this 
orientation. Besides, the angular selectivity of point con- 
tact measurement is not high enough. Since tunneling 



FIG. 9: Temperature dependences of the broadening F 
(closed circles) and pair-breaking 7 (open circles) parame- 
ters for LuNi2B2C-Ag point contact in the ab plane (a) and 
in the c direction (b) in the two-gap approximation. Note 
that the parameters are considerably higher for the small gaps 
(squares). 



component in the point contacts is rather weak (the Z 
coefficient is not high) and their diameters are close to 
the coherence length, it is quite likely that the conversion 
of the electrons into pairs in the vicinity of the contact 
occurs within a spherical geometry. According to Maki's 
modelS^ the angle between the largest and the smallest 
gaps in the ab plane is 45°. As a result, the spectra taken 
in these directions can hardly have big distinctions. The 
most obvious features discovered by us are the different 
critical temperatures of the small gap in the ab plane and 
in the c direction obtained in the BCS extrapolation (Fig. 
0. 

With the evidence available, it is hardly possible to 
conclude unambiguously whether this phenomenon is re- 
lated to the gap anisotropy or any other (e.g., technolog- 
ical) factor. In the ab plane the point contact was formed 
on a freshly cleaved surface. In the c direction the sur- 
face was cleaned chemically and experienced no mechan- 
ical stress. It is known that elastic scatterers first cause 
gap isotropization, then reduce the gap and suppress the 
critical temperature. Although in our experiments the 
critical temperatures of the point contacts coincide with 
Tc of the bulk, the partial Tc of the small gap ('--' lOK Fig. 
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FIG. 10: Temperature dependences of the rms deviation 
characterizing shape discrepancies between experimental and 
theoretical curves (calculation by the equations of Ref.i2l| 
(open circles) and in the BTK model Ref.jgO (solid circles)) 
for LuNi2B2C-Ag point contacts in the ab plane (a) and in 
the c direction (b) in the two-gap approximation. 



I?)) and its relative contribution (coefficient Ki Fig. ISJ to 
the general conductivity decreases in the ah plane, which 
indicates isotropization of the general gap in this plane. 
To clear up possible reasons for this phenomenon, fur- 
ther measurements on chemically cleaned natural growth 
faces of single crystals in the a direction are needed. 

According to Pokrovsky's theorem22i2S in the case of 
anisotropic one-gap superconductivity, the temperature 
dependence of the gap is independent of crystallographic 
directions. In our case the temperature dependences of 
the large and the small gaps and the BCS extrapola- 
tion of the critical temperatures are quite different. This 
prompts the conclusion that this behavior can not be 
caused by the contributions to conductivity made by a 
single anisotropic superconducting gap in different crys- 
tallographic directions. On the contrary, we observe the 
indications of two different superconducting gaps from 
different energy bands. 

Noteworthy also is another point, namely, the limits of 
applicability of the two-gap approach to describing ex- 
perimental data. In the strict sense, the method requires 
that the broadening (pair breaking) parameters should 
be zero for both gaps. This condition is met at T > Q K 



for ah plane and T > 8 ii' in the c direction (see Fig. 
At lower temperatures the gaps do not have strictly spec- 
ified values as they are distributed over a certain range 
of energies simulated by nonzero F or 7. In the ah plane 
the largest gap stands out as an individual line above 2 
K , while in the c direction both the gaps start broaden- 
ing at nearly the same temperature. As the temperature 
lowers, the distribution extends up to overlapping of en- 
ergies. This brings more uncertainly into the description 
of curves in the two-band approximation (Fig. I10|l and 
sends us in search of a new method. A possible tech- 
nique of a more adequate description of our experimental 
results near T = l.h K is illustrated in^ (Figs. 12-14). 



VI. CONCLUSIONS 

This study demonstrates that the two-gap approxima- 
tion is more appropriate for describing the superconduc- 
tivity of LuNi2B2C in a wide interval of temperatures. 
The values and the temperature dependences of the large 
and the small gaps have been estimated for ah plane and 
the c direction. It is found that in the BCS extrapolation 
of the critical temperature for the small gap is Tc = 10 
K in the ah plane and Tc — 14.5 K in the c direction. 
The absolute values of the gaps are Ag* = 2.16 meV, 
Aq = 1.94 meV . For the large gaps Tc coincides with 
j^buik. 7;(A2) 16.8 K and their absolute values are 
very close in both orientations, being about 3 meV. It is 
found that in the c direction the contributions to conduc- 
tivity from the small and the large gaps are nearly equal 
up to 10^11 K. In the ab plane the contribution of the 
small gap is considerably weaker and decreases rapidly 
as the temperature is growing. 
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APPENDIX A 

The iteration method commonly used to fit theoretical 
and experimental curves is quite good when broadening 
is small. However, as the broadening F becomes compa- 
rable with the gap (or 7 > 0.3), this introduces an appre- 
ciable uncertainly into the results. The iteration method 
implies that during fitting the parameters change cycli- 
cally as the step is gradually decreased. It is expected 
that the error reduces with each step. Our experiment 
shows that with high F (or 7) this procedure brings us 
finally to a local minimum point which is dependent on 
starting A, F (7) and Z values and cannot correspond 
to the global minimum. Moreover, this method is not 
valid for the two-gap calculation. In our fitting proce- 
dure we therefore used the technique of coordinate de- 
scent with a postponed solutionis. First, we specified an 
interval in which A is searched for at a given tempera- 
ture. The interval was then subdivided into equidistant 
parts Ai, A2, An. Then F (7) and Z were fitted for 
each Ai. The method used in^^ does not imply that 
the error should decrease immediately after each F (or 
7) step. The errors are compared after each F (or 7) 
step only when the Z-fitting is completed (this accounts 
for the term "postponed solution"). In the programming 
language, the calculation by this method reduces to em- 
bedded cycles. In the one-gap model, the Z-fitting is an 
inner cycle and the F (or 7) fitting is an outer cycle. In 
the two-band calculation, the Z-fitting is the innermost 
cycle, the fittings in Fi and F2 become outer and, finally, 
the A2-fitting is the outermost cycle. This method pro- 



vides an unambiguous solution independed of the starting 
parameters. In the process of calculation the intensities 
of the theoretical and experimental curves were equalized 
and aligned over the ordinate before each step of calcu- 
lation of the average rms deviation i^(Ai). To do this, 
the y-coordinates of the points in the theoretical curve 
were multiplied by the scale factor S (the procedure of 
making the amplitudes of the theoretical and experimen- 
tal curves equal). The curve was then shifted along the 
y-axis by amount B. The values of S and B were found 
from the condition of the minimum ^(Aj). The standard 
algorithm for determination of these coefficients known 
as the least-square method is considered, for instance, 
inSli. As a result, for each A; we could obtain F^ (7^) 
and Z at which the difference between the shapes of the 
theoretical and experimental curves characterized by the 
rms deviation _F(Ai) was the smallest one. The kind of 
calculation for different temperatures is illustrated iniS 
(Fig. 3). In the two-gap case, the i^(Ai)-curves are more 
fiattened, and the curve has no minimum at a tempera- 
ture lower than in the one-gap case. However, before this 
happens, we are able, as a rule, to find quite accurately 
the scale-factor S which is independent of temperature 
but is strongly dependent on A (— , Fig. 4). The value 
of S can vary from contact to contact but it is invari- 
ant for a particular contact. To put it more accurately, 
at the values of F ~ 0, 5' is practically invariable for 
the minima in the curves i^(A) (see Fig. 9 inl?). By 
choosing A-values which correspond to the determined 
magnitudes of S*, we can plot quite accurately the curve 
A(T) in the high-temperature region. 
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